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Abstract 

The problem of wave disturbance propagation in rarefied gas in gravity field is 
explored. The system of hydrodynamic-type equations for a stratified gas in gravity 
field is derived from BGK equation by method of piecewise continuous partition 
function. The obtained system of the equations generalizes the Navier-Stokes at 
arbitrary density (Knudsen numbers). The verification of the model is made for a 
limiting case of a homogeneous medium. Results are in the good agreement with 
experiment and former theories at arbitrary Knudsen numbers. 

Introduction 

There is a significant number of problems of gas dynamics at which it is necessary to use 
the mathematical apparatus beyond the limits of traditional hydrodynamics of Navier - 
Stokes. The hydrodynamics is valid under the condition for the Knudsen number Kn = 
l/L <C 1, where / is a mean free path, and L is a characteristic scale of inhomogeneity 
of a problem under consideration. The first work, in which wave perturbations of a gas 
were investigated from the point of view of more general kinetic approach, perhaps, is the 
paper of Wang of Chang and Uhienbeck pQ. Most consistently these ideas are formulated 
in the work of Foch and Ford j2j- Such general theory could be based on some kinetic 
approach, i.e. Boltzman equation. 

Numerous researches on a sound propagation in a homogeneous gas at arbitrary Knud- 
sen numbers were made j3] - ^2] • The investigations have shown, that at arbitrary Knud- 
sen numbers the behaviour of a wave differs considerably from ones predicted on a basis 
of hydrodynamical equations of Navier - Stokes. These researches have revealed two es- 
sential features: first, propagating perturbations keep wave properties at larger values of 
Kn, than it could be assumed on the basis of a hydrodynamical description. Secondly, at 
Kn > 1 such concepts as a wave vector and frequency of a wave become ill-determined. 

The case, when Knudsen number Kn is non-uniform in space or in time is more 
difficult for investigation and hence need more simplifications in kinetic equations or their 
model analogues. A constructions of such approaches for analytical solutions based on 
kinetic equation Bhatnagar - Gross - Krook (BGK) of Gross- Jackson [3] in a case of 
exponentially stratified gas were considered at [2B1 EH] 

In this paper we would develop and generalize the method of a piecewise continuous 
partition functions [211121] to take into account the complete set of nonlinearities (?). We 
consider the example of wave perturbations theory for a gas stratified in gravity field so 
that the Knudsen number exponentially depends on the (vertical) coordinate. 
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1 Piecewise continuous partition function method 



The kinetic equation with the model integral of collisions in BGK form looks like: 

df ^df df 

here / is the distribution function of a gas, t is time, v is velocity of a particle of a gas, r 
is coordinate, 

^ * n ( (v-Uf 
fi[r,v,t)) = -^y^exp — 

is the local-equilibrium distribution function, H = kT/mg is the so-called height of a 
homogeneous atmosphere - a parameter of the gas stratification, v T = \j2hT /m is the 
average thermal speed of movement of particles of gas, v = u exp(—z/H) is the effective 
frequency of collisions between particles of gas at height z. It is supposed, that density 
of gas n, its average speed U = (u x , u y ,u z ) and temperature T are functions of time and 
coordinates. 

Following the idea of the method of piecewise continuous distribution functions let's 
search for the solution / of the equations (j2J) as combinations of two locally equilibrium 
distribution functions, each of which gives the contribution in its own area of velocities 
space: 



/ m \ 3 / 2 



f{t,f,V)= I 




ffl \ 3/2 



(2) 



here n , U ± ,T ± are parameters of locally equilibrium distributions functions. Geometry 
of break, that is the area, in which various functions operate, is determined by geometry 
of a problem. 

Thus, a set of the parameters determining a state of the perturbed gas is increased 
twice. The increase of the number of parameters of distribution function (J2J) results in 
that the distribution function generally differs from a local-equilibrium one and describes 
deviations from hydrodynamical regime. In the range of small Knudsen numbers I « L 
we have n + = n~,U + = U~,T + = T~ and distribution function (0) tends to local- 
equilibrium one, reproducing exactly the hydrodynamics of Navier-Stokes. In the range 
of big Knudsen numbers the formula (J2J) gives solutions of collisionless problems. Similar 
ideas have resulted successfully in a series of problems. For example, in papers [IHj - [2H1 a 
method of piecewise continuous partition function was demonstrated for the description of 
flat and cylindrical (neutral and plasma) Kuette flows JH] - [2D] • Thus for a flat problem 
the surface of break in the velocity space was determined by a natural condition V z — 0, 
and in a cylindrical case V r — 0, where V z and V r are, accordingly, vertical and radial 
components of velocity of particles. Similar problem was solved by perturbations caused 
by pulse movement of plane [201 121]- Solving a problem of of a shock wave structure (201 
1221 122] the solution was represented as a combination of two locally equilibrium functions, 
one of which determines the solution before front of a wave, and another - after. In the 
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problem of condensation/evaporation of drops of a given size |24[ 125] a surface break was 
determined by so-called "cone of influence" , thus all particles were divided into two types: 
flying " from a drop" and flying " not from a drop" . 

The similar approach was developed for a description of a nonlinear sound in stratified 

gas num. 

The idea of a method of two-fold distribution functions given by (j2J) is realized as 
follows. Let's multiply equation BGK (0) on a set of linearly independent functions. In 
the one-dimensional case U = (0, 0, U z ) the following set is used: 



(3) 



fi = rn , 


ip 4 


= m(y z - 


u z )\ 


cp 2 = m V z , 


¥5 




)\V -u\ 2 


<^3 = \m\V-U\\ 




\ m{ y z - 


- u z f . 


scalar product: 








< fn, f > = < <Pn > ; 


'J 


dvip n (t,z,V)f(t. 


,z,V) . 



(4) 



< <pt >= m < 1 >= p(t, z) , <cp 4 >=m< (V z - U z f >= P zz , 

< (p 2 >=m < V z >= pU z , < v^s >= < (V z - U Z )\V-U\ 2 >= q z , ( 5 ) 

< (p 3 >= l -m < \V - U\ 2 >= e , < <p 6 >= -m < (V z - U z f >= q z . 

Here p is density, pU z is specific momentum, e is internal energy per unit mass of the 
gas, P zz = P + n zz is the diagonal component of the strain tensor (P is pressure, tt zz 
is component of strain tensor), q z is a vertical component of a heat flow, q z is the new 
parameter having meaning of a heat flow. 

Multiplying equation (0) on eigen functions © we obtain the system of differential 
equations: 

d d 1 d 

7 rU z + U z —U z + - — (P + U zz )+g = 0, 

at oz p oz 

—e + U z —e +{e + P + U zz )—U z + —q z = , 

(Jv CJZ CJZ CJZ f n\ 

-{P + U zz ) + U Z —(P + U zz ) + 3(P + U zz )—U z + 2—q z = -u(IL„ + P-p9), 

at oz oz oz 

d d d 1 d d 

—q z + U z —q z + 2{q z + q z )^-U z - -(e + P + n,,) — (P + U zz ) + —J x = -vq z , 

at oz oz p oz oz 

I*, + U z j-q z + 4q z §-U z - |-(P + U zz )j-(P + H„) + I J, = -uq z , 

where J x =< (V z - U z f{V - Uf > , J 2 =< (V z - U z f > . (7) 
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The obtained system © of the equations according to the derivation scheme is valid at 
all frequencies of collisions and within the limits of high frequencies should transform to 
the hydrodynamic equations. It is a system of hydrodynamical type and generalizes the 
classical equations of a viscous fluid on any density, down to a free-molecule flow. However, 
the system (JOJ) is not closed yet. It is necessary to add equations of state P = P(p,T) 
and e = e(p,T). Except for that it is necessary to present values of two integrals J\ and 
J2 as functions of thermodynamic parameters of the system © . 

Let's evaluate integrals (J7J) directly, plugging the function ©. We estimate the func- 

tions — as small, that corresponds to small Mach numbers M = max\v\/vT- Values of 
V T 

integrals J\ and J2 within the specified approximation looks as 

■h = 4<" +l/ T +4 + n-Vf 4 ) + -^[n+V'+ 3 (f/+ -U)- n-V^iU- - U)\ , 
J 2 = _( n +i/+ 4 + n~Vr 4 ) + -^[n + V+ 3 (U + -U)- n~Vf 3 (U- - U). 

16 \/1t 



Let's express parameters of the two-fold distribution function (J2J) through the thermody- 
namic ones and substitute the result into the expression (JHJ). To solve the specified problem 



we shall use a method of perturbations with the small parameter ma ^ - ; expressing 



U + = vf (Pzz - fe) + i(5Ue + 3q) 
U- = ^J2(P ZZ ~ h) + hi^Ue + 3q), 



p- = P+ U P » - ¥) - l^(f) 3/2 ^ ( 9 ) 



3p 12 Y pe 
^ 4- — . 

3p~12Vpe^ zz 3"^ 1 10 e 

Plugging the values of © into (JHJ) one obtain the values of Ji j2 in the first order : 



10e 2 61e/ 2\ T 2 e 2 13 e / 2\ 

2 Limiting case of gas oscillations at high frequencies 
of collisions (small Knudsen numbers). 

Let us consider a system in the hydrodynamical limit [v — > 00). It follows from the last 
three equations of the system (jUJ) that the orders of values relate as max{U zz ,q z ,q z } ~ 
u^ 1 max{p, U z , e, P}. Next assume v~ x = in the zero order by the parameter v~ x . One 
hence have J\ zz = 0, q z = 0, q z = at the l.h.s. and at the r.h.s. of forth equation of the 
system P = p9 = 2e/3. Substituting mentioned limits in the first three equations of the 
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system © we obtain a system of Euler equations of a liquid in gravity field: 

§i u - +u -ii u - + hi? + » =0 

d TT d 5 d TT 
— e + f/ z — e + -e— t/ z = . 

at (72; 3 dz 

The functions {Il zz ,g 2 ,g z } ~ u~ 1 {p,U z ,e, P} belong to the next order of the parameter 
u^ 1 . Then from the last three equations of the system, © taking into account the equation 
of state P = p9, one obtains following relations 

8 d TT 10 d fe\ _ 2 d f e\ 

n zz = -7r-r^e—U z , q z = —p—r^^- - , Qz = rrTT - ■ ( 12 ) 



9v(z) dz ' 9v(z) dz \p J ' 3v(z) dz \p y 

Further substituting (fT2*|) in the first three equations of the system (jOJ) we obtain 

d TT TT d TT Id 8 d (e d \ 

df* + Uz d-z U > + 3-pd-z e + 9 -Ypdz \udf z ) = ° ( 13 ) 

d TT d 5 d TT 10 d (e d e\ 8e / d r _ , 
— e + C4— e + -e—U z - — — ( - — - ) - -- ( —U z ) = . 



dt dz 3 dz 9 dz \v dz p J 9 u \ dz 

System ()13)1 is the system of equation of a non-ideal liquid, to compare it with the Navier 
- Stokes equations, we continue the evaluation of viscosity factor and coefficient of heat 
conductivity. Expressions for strain tensor and heat flow tensor in one-dimensional hy- 
drodynamics take a form 

4 9 n 2 d n 

Kzz = V-T^Uz , Qz = ^TTP , 

3 'dz 3 dz 

where rj is the viscosity factor, and k is the coefficient of heat conductivity. Comparing 
mentioned expressions with corresponding items in equations f!13|) we obtain 

n kT 5n kT , 

v = -—> K = - 2 —' ^ 

that coincides with the well known relations, given, for example, in [3U]. Finding the 
Prandtl number, taking into account of molecular thermal capacity of the ideal gas under 
constant pressure C p = 5/2, we obtain 

Pr = ^ = l, 



that do not coincide with the Prandtl number of ideal gas (Pr^ = 2/3). The wrong 
Prandtl number is the main disadvantage of BGK model, that, however, can be removed 
by changing to the more exact models of collision integral, for example, of Gross- Jackson 
0- 



5 



3 Linearized system of the equations. Dispersion re- 
lation. 



(15) 



For a closure of the system we use the equation of state of ideal gas. Linearized system 
of the equations (|TT|) is given by 

d d 

|^ + 5^ + T + n -) = ' 

^+1^(2^ + 3^=0, 

d 1 d 

^n 22 + -V T —{AU Z - 3q z + 9q z ) = -uU zz , 
at 3 oz 

d X d 

Ft q * + 36^ (30T + 3in " ) = -" 9 " • 
| ? , + _Lv4(6T+7n«) = -^. 

For convenience we would introduce new notations Hi for hydrodynamical variables: n\ = 
p, n 2 = U z , n 3 = T, 714 = U zz , n 5 = q z , n 6 = q z . The solution of system (fTo]) we search as 

rii = diexp(—iwt + ik z z), (16) 

where w is frequency of a wave, k z - the vertical component of a wave vector. 

Substituting (fTo^) in (|T5*j) .one obtains a system of the homogeneous algebraic equations 
with constant coefficients which solution exists if 

18 ~ fi /3 2 39 48 \~ 4 / o 24 2 23 58\ ~ 2 , 2 
— k b + -r 2 - — - — ir k i + -ir s - —r 2 + —ir + — k +ir —l — 3ir + 3r = 

IzO /V o ID / 

(17) 

Here the dimensionless wave number k = kCo/w and the Reynolds number r = vjw are 
introduced, where Co = yf^r - sound speed in Euler's approximation. The Reynold's 
number r and the Knudsen number are obviously linked: 

A w Vr /6 1 

Kn 



Xb v 2ttCo V 5 27rr 
Let k = (3 + ia, then 

Hi = aiexp(—tw{t — — z))exp(— w— z) 

Cq L'O 

and the real part /3 = Cq/C, a - the factor of attenuation. 

4 The joint account of three modes. 

The basic Fourier component solution of the system (|THj) we shall search as a superposition 
of three plane waves 

rii = A]exp(—iwt + ik\z) + A 2 exp(—iwt + ik 2 z) + A^exp(—iwt + ik%z), (18) 
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where kj, j = 1,2,3, are solutions of the dispersion equation (JTTj) correspondent to the 
modes. 

Substituting (fTHj) into the linearized system, we express A J 2 , A 3 , A{, A 3 5 , A\ through 
A{ = A j . For A 3 2 , A J 3 we have: 



Ai wAj 



2 ' 

j A 1 (-31 V T 2 k 2 + 2Uw 2 - 2Avw + Q2w 2 ) 

3 ~ V%k) + 36iw 2 - 36uw 



(19) 



To determine the coefficients A 1 , A 2 , A 3 we should choose boundary conditions. We con- 
sider a problem in half-space and the reflection of molecules from a plane as a diffuse one 
[7|. The boundary condition for the distribution function looks as 



f( z = 0,V,t) = -^exp{- {V ^ J -} bjV z >0. 



-iwt\2 



Here Uq stands for an amplitude of the hydrodynamic velocity oscillations. For ^ < 1 
we have: 

<p(z = 0, V, t) = LzJ^L „ 2^e~ M by V, > 0. 
For hydrodynamical variables on the bound we obtain 

p ( z = o, t) =< <p( z = o, v, t) >= [ dv^{z = o, v, ty- v2 i y T = ^= e ~ iwt 



IK ' V rp 

U z (z = 0, t) =< ^f<p(z = 0, V, t) >= ^e~ M 
T(z = 0, t) =< ^<p{z = 0, V, t) >= -^e"*" 

Vrr WIT 



T 



(20) 

Substituting the values of (fTH|) into (fTHj) and comparing right-hand sides of expression 
(j!8)l and (J20J) we obtain the system of equations in variables AK Solving given system of 
equations we obtain variables A 1 , A 2 , A 3 . 

In experiment acoustic pressure perturbation amplitude is measured. Appropriate 
combination of the basic variables for the pressure it is given by the formula 

P(z, t)e iwt = p' + T'= {A 1 + A\)e ikiz + (A 2 + A 2 3 )e ik2Z + {A 3 + A\)e ik3Z (21) 

The real part of this expression relates to experiment. In fig. la) the real part of this 
expression is represented at r=0.2, where Z = zw/Cq - dimensionless coordinate. The 
attenuation factor a is determined as a slope ratio of the diagram of the logarithm of 
amplitude of pressure depending upon distance between oscillator and the receiver. It is 
illustrated in fig. 1 b). 
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Fig. 1. Finding of attenuation factor 

5 Comparison with with experimant and results of 
other evaluations. 

In figures 2,3 a comparison of theoretical results of the sound propagation parameters 
with experimental data |U E] is made. 

The dispersion relation (j!7j) represents the binary cubic equation with variable coeffi- 
cients. The exact analytical solution by the formula Cardano is very huge and therefore we 
do not show it in this paper. At r — > (free molecule flow) we start from the propagation 
velocity by the formula 

C /C = 0.54 + 0.15r 2 + 0(r 4 ) 
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Fig.2. The propagation velocity. Thick line - Navier-Stokes, thin line - this paper, circles 
- measurements in Argon (Greenspan, Meyer-Sessler) . 



The attenuation factor is determined graphically as shown in the fig. 1. Therefore we 
cannot introduce analytical expression. 




0.01 0.1 r 10 



Fig.3. Attenuation in Argon, thick line - Navier-Stokes. thin line - this 
paper (sound- wave), dotted line - the joint account of three modes 

Results for phase speed give the good consistency with the experiments. As we see, 
the account of three modes allows us to enter further the area of intermediate Knudsen 
numbers. 

In figures 4, 5 a comparison of our results of numerical calculation of dimensionless 
sound speed and attenuation factor depending on r is carried out with the results of the 
other authors. 
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Fig. 4. A comparison of BGK model results with experimental data. The propagation 
velocity. 1 - Present work, 2 - Loyalka and Cheng JI],3- Buckner and Ferziger [7] ,4- 
Navier-Stokes theory, 5- Sirovich and Thurber j^J, 6-experimental data of Greenspan, 7- 

experimental data of Meyer and Sessler 



0.4 1 




Fig. 5. A comparison of BGK model results with experimental data. The attenuation 
factor. 1 - Present work, 2 - Loyalka and Cheng [TT], 3- Buckner and Ferziger 4- 
Navier-Stokes theory, 5- Sirovich and Thurber j^J, 6-experimental data of Greenspan, 7- 

experimental data of Meyer and Sessler 

Lacks of the BGK model used in this article is that it gives correct value of viscosity 
factor, but wrong value of coefficient of heat conductivity. To the superior models of the 
Gross - Jackson ^3J this lack can be eliminated by transition. 

At the solution of Boltzman equation the method of the Gross - Jackson revealed 
sudden disappearance of discrete modes at some values r c (0, [7j, jUj), and with increase 
of number of the moments r c decreased. 

For example, Buckner and Ferziger in the paper [7j have shown, that for r > 1 the 
solution is determined mainly by the discrete sound mode and the dispersion relation 
may be used in calculating the sound parameters. For r < 1, the continuous modes are 
important. The solution remains "wavelike", but it is no longer a classical plane wave. In 
fact, the sound parameters are depend on the position of the receiver. 

Below r c the solution is represented as superposition of a continuous spectrum of 
eigen functions, therefore the classical understanding of a sound should be changed. The 
concept of a dispersion relation is not applicable more. 

6 Conclusion 

The attenuation of sound at big Knudsen numbers is not " damping" (due to intermolecular 
collisions), but rather "phase mixing" (due to molekules which left the oscillator at different 
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phases arriving at the receiver at the same time) 

The attenuation factor at big Knudsen numbers Kn > 1 is modelled by the account of 
effects of a relaxation in integral of collisions. The model of the Gross - Jackson at given 
N limits an opportunity of the account external times of a relaxation (fast attenuation) 
as essentially bases on a condition: 

A» = Ajv+i, i>N + l 

Supreme times of a relaxation are assumed identical. It means, that the inclusion of 
the supreme eigen functions \ii i > N + 1 is necessary, that would allow to move in 
the range of higher Knudsen numbers. 

In piecewise continuous partition function method the number of waves is twice more, 
but restrictions on attenuation factor remain. 
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